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Limits on v u 



and for Am 



i/ e and v 
M beams. 
> 1000 eV 2 , 



—* v e oscillations are extracted using the NuTeV detector with sign- 
In Tip. mode, for the case of sin 2a = 1, Am 2 > 2.6 eV 2 is excluded, 
sin 2 2a > 1.1 x 10" 3 . The NuTeV data exclude the high Am 2 end of 
Dp — > F e oscillations parameters favored by the LSND experiment without the need to assume that 
the oscillation parameters for v and V are the same. We present the most stringent experimental 
limits for — > v e (V e ) oscillations in the large Am 2 region. 



PACS numbers: 14.60.Pq, 13.15.+g ; UR-1640 to be published in Phys. Rev. Lett. 



Experimental evidence for oscillations among the three 
neutrino generations has been recently reported. For 
two-generation mixing, the probability that a neutrino 
created as type v\ oscillates to type is: 



Piv\ — ► v-i) = sin 2 2a sin 2 



1.27Am 2 L 



(1) 



where Am 2 is the mass squared difference between the 
mass eigenstates in eV 2 , a is the mixing angle, E v is the 
incoming neutrino energy in GeV, and L is the distance 
between the points of creation and detection in km. 

Data from the Super-Kamiokande atmospheric neu- 
trino experiment Q have been interpreted as evidence for 

— > v T oscillations with sin 2 2a > 0.88 and 1.6 x 10~ 3 < 
Am 2 < 4 x 1CT 3 eV 2 . The LSND experiment has re- 
ported jf|] a signal consistent with — > v e oscillations 
with sin 2 2a « 10~ 2 and Am 2 ~ 1 eV 2 . The solar 
neutrino experiments, and most recently SNO have 
reported evidence for oscillations of v e — > (V m ,j/ t ) with 
Am 2 < 10 _3 eV 2 . Within a three-generation mixing sce- 
nario and under the assumption that the Am 2 values for 
v and V are the same, it is not possible to simultaneously 
accommodate the Super-Kamiokande, LSND, and SNO 
results. Therefore, experimental searches for oscillations 
with both v and v beams are of interest. In this letter, 
we report on a search for oscillations in both the v„ — > v e 



High-purity v and V beams are provided by the new 
Sign-Selected Quadrupole Train (SSQT) beamline at the 
Fermilab Tevatron during the 1996-1997 fixed target run. 
Hadrons are produced when the 800 GeV primary proton 
beam interacts in a BeO target located 1436 m upstream 
of the neutrino detector. Sign-selected secondary parti- 
cles of specified charge (mean energy of about 250 GeV) 
are directed in a 221 m beamline towards a 320 m decay 
region, while oppositely charged (and neutral) mesons 
are stopped in beam dumps. Two-body decays of the 
focused pions yield (77^) with a mean energy of «75 
GeV. Two-body decays of the focused kaons yield (V^) 
with a mean energy of w200 GeV. Muons are stopped in 
a 915 m steel/earth shield. 



The energy and spatial distributions of 



cc 



and 
beam. 



channels using a new sign-selected neutrino 



events in the detector provide a determination of the 
flux of pions and kaons in the decay channel (used in 
the determination of the predicted v e and v e fluxes). 
For running mode, the predicted energy spectra for 
Vu, Vfj,, and (z^ e +77 e ) CC events are shown in Figure 
|l|(a). The corresponding spectra for V ^ running mode are 
shown in Figure 0(b). The (7/ M )beam contains a 1.7% 
(1.6%) ^ e 's(lVs) 93% and 70% of which are produced 

from K — > 7r°e ± vl, for v and v modes, respectively. 
The proton beam is incident on the production target at 
an angle such that forward neutral kaons do not point at 
the detector. This greatly reduces the electron neutrino 
flux from neutral kaon decays (which is more difficult to 



1 



model) . The error in the predicted electron neutrino flux 
is reduced from 4.1% (in CCFR j§) to 1.4% (NuTeV). 

The NuTeV detector |] is an upgrade of the CCFR 
detector [F|. It consists of an 18 m long, 690 ton to- 
tal absorption target-calorimeter with a mean density 
of 4.2 g/cm 3 . Muon energy is measured by a 10 m 
long iron toroidal spectrometer. The target consists of 
168 steel plates, each 3 m x 3 m x 5.15 cm, instrumented 
with liquid scintillation counters placed every two steel 
plates and drift chambers every four plates. The sepa- 
ration between consecutive scintillation counters corre- 
sponds to six radiation lengths. The energy resolution 
of the target calorimeter is AE h /E h w 0.85/-\/E*"(GeV), 
and AE e /E e w 0.50/\/5e(GeV) for hadrons and elec- 
trons, respectively. The muon momentum resolution is 



0.11. The NuTeV detector is calibrated con- 



tinuously every accelerator cycle (once a minute) with 
beams of electrons, muons, and hadrons during the slow 
spill part of the cycle. 
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FIG. 1. (a) The predicted visible energy spectra for , 77 M , 
and (v e ,V f ,) CC-events in running mode, (b) The corre- 
sponding spectra for 77^ running mode. The predictions come 
from a beam simulation tuned to agree with the observed 
number of or 77 M CC events in each running mode. 

While the neutrinos arrived in gates a few msec wide, 
the calibration beam arrived in a different gate 20 sec 
long, followed by an off-spill cosmic ray gate for back- 
ground measurement. These continuous test beam cal- 
ibrations yield a reduction in the hadron energy scale 
error from 1% (in CCFR §) to 0.43% (in NuTeV §). 

The event sample used in this analysis is similar to 
that used in the recent precise NuTeV measurement of 
the electroweak mixing angle J7J with additional fiducial 
cuts, and E ca \ > 30 GeV. The data sample consists of 
1.5 x 10 6 v events and 0.35 x 10 6 V events with a mean 
visible energy in the calorimeter (i? C ai) of 74 GeV and 



56 GeV, respectively. The observed neutrino events are 
separated into CC and NC candidates. Both CC and 
NC interactions initiate a cascade of hadrons in the tar- 
get that is registered in both the scintillation counters 
and drift chambers. Muon neutrino CC events are dis- 
tinguished by the presence of a final state muon, which 
typically penetrates well beyond the hadronic shower and 
deposits energy in a large number of consecutive scintil- 
lation counters. NC events usually have no final state 
muon and deposit energy over a range of counters typ- 
ical of a hadronic shower (about ten counters ~ 1 m of 
steel). For each event, the length (L) is defined as the 
number of scintillation counters between the interaction 
vertex and the last counter consistent with at least single 
muon energy deposition. A pure sample of v^N — > fj,~X 
charged current events is obtained from a 'long' sam- 
ple with L > 29 for v running mode (L > 28 for V). The 
"short' event sample consists of events with L < 28 for 
v running mode (L < 27 for v). Events with a 'short' 
length arc primarily NC induced and originate from: 



1. V»,eN 



i/» e NC events; (« 65%); 



2. v^N — > n X\ short CC events with muons 
which range out or exit the side of the calorime- 
ter (25% for v, 15% for V); 

3. v e N -> eX v e CC events (10% for u, 15% for 77); 

4. (j,N — > fiX; steep cosmic ray interactions (2% for v 
and 9% for 77). 

The electron produced in a v e CC event (source 3) de- 
posits energy in a few counters immediately downstream 
of the interaction vertex; this changes the longitudinal 
energy deposition profile of the shower. The energy pro- 
file is characterized by the ratio of the sum of the energy 
deposited in the first three scintillation counters to the 
total visible energy in the calorimeter E ca \: 



Ei + Ei 



E, 



E, 



(2) 



where E t is the energy deposited in the i th scintillation 
counter downstream of the interaction vertex, and E c& \ 
is the sum of the energy in the 20 scintillation counters 
downstream (plus 1 upstream) of the vertex. We simi- 
larly define 772 to be the ratio of the sum of the energy 
deposited in the first two scintillation counters to E ca \. 

Although the electromagnetic shower component is 
typically much shorter than a hadron shower, NC and 
v e CC interactions cannot be separated on an event-by- 
event basis. However, the difference in the shower pro- 
files can be used to perform a statistical extraction of the 
number of v e CC events in the "'short' sample using the 
773 distribution. 

We first assume that (for the same final state hadron 
energy) hadron showers produced in NC and CC inter- 
actions are the same. Any difference in the 773 (or 772) 



2 



distributions of 'long' and 'short' events is attributed to 
the presence of v e CC interactions in the 'short' sam- 
ple. To compare directly the 'long' and 'short' events, a 
muon track from the data is added to the 'short' events to 
compensate for the absence of a muon in NC events. The 
fraction / of CC (source 2) events with a low energy 
muon contained in the 'short' sample which now have two 
muon tracks is estimated from a detailed Monte Carlo of 
the experiment. A simulated sample of such events is 
obtained by choosing 'long' events with the appropriate 
energy distribution from the data to which a second short 
muon track is added in software. The length of the short 
track and the angular distribution are obtained from a 
Monte Carlo of CC events. 

A sample of v & CC interactions in our detector is simu- 
lated by adding a GEANT (To) generated electromagnetic 
shower of the appropriate electron energy to the calorime- 
try counter energies in events in the 'long' data sample. 
There is good agreement between the GEANT simula- 
tion of electrons and the test beam data. The energy 
distribution of the electron neutrinos and the fractional 
energy transfer y (in each event) are generated using a 
detailed Monte Carlo simulation of the experiment. The 
CC cross section model is tuned to agree with the mea- 
sured CCFR differential cross sections Q for events. 
Since the hadron showers in the 'long' sample already 
have a muon track, this sample ( v e CC + fi) can be com- 
pared directly with the 'short' and 'long' samples. 

The 'long' and 'short' 773 distributions are further cor- 
rected by subtracting the contamination due to cosmic 
ray events. The cosmic ray component (source 4) , which 
is only important for very low energy bins, is well mea- 
sured using the off-spill cosmic ray data. Additionally, 
the 773 distribution of short CC events (source 2), nor- 
malized to the predicted fraction /, is subtracted from 
the 'short' event sample. 

To extract the number of v & CC events in each E ca \ bin, 
we fit the corrected shape of the observed 773 distribution 
for the 'short' sample to a combination of CC and 
v e CC distributions with appropriate muon additions: 



[short +fi\=a [f^CC] + (3 [v c CC + Li] 



(3) 



At this point we improve || over the previous CCFR 
analysis [Q by correcting for additional effects. First, 
the hadron shower in CC events also includes the con- 
tribution of photons that are radiated from the muon 
during the CC scattering process. These photons are not 
present in hadron showers of NC events. A correction is 
applied by using the PYTHIA Monte Carlo Q to gen- 
erate the spectrum of photons radiated by the muon in 
CC events. The parameters in PYTHIA that govern the 
emission of photons are tuned to yield agreement with 
the radiative corrections formalism of deRujula |Q. Sec- 
ond, the procedures of adding a muon in software to the 
'short' sample (to model a 'long' event) and of modeling 



electron neutrino events by adding GEANT-generated 
electromagnetic showers to hadron showers from 'long' 
events are corrected for imperfect modeling using a full 
LEPTO /GEANT/GHEISHA simulation of the experi- 
ment and analysis procedures. In the GEANT simula- 
tion of neutrino events, the Lund Model is used to gen- 
erate the initial particle composition of hadron showers. 
The entire experimental procedure is simulated with the 
beam Monte Carlo (i/ e ,V e ) flux as input. Modeling cor- 
rections are extracted for each E ca \ bin from the small 
difference between the extracted (v e ,V e ) flux using simu- 
lated data and the flux extracted using perfect modeling 
in the Monte Carlo. 

The absolute flux of v^s is taken as the average of the 
results from analyses done using the 773 and 772 variables 
(the statistical error from the 773 analysis is used) . An 
additional systematic error (2.3% in neutrino mode and 
0.6% in antineutrino mode) is included to account for 
this difference. This systematic error for antineutrinos is 
smaller because the final state positron carries a much 
larger fraction of the energy in V e CC events. 
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FIG. 2. The ratio of the detected over predicted num- 
bers of (vepe) events versus visible energy minus 1. The 
curves correspond to the predictions for oscillations with 
sin 2 2a = 0.01, and Am 2 of 100 and 1000 eV 2 . The solid 
line is the 90% confidence upper limit for Am 2 =1000 eV 2 . 
The shaded area corresponds to the systematic error band. 



For the oscillation search we measure the absolute flux 
of v e J s at the detector and compare it to the flux pre- 
dicted by a detailed beamline simulation 0(9). In order 
to extract limits on oscillations, the data are fitted by 
forming a x 2 which incorporates the Monte Carlo gener- 
ated effect of oscillations and terms with coefficients ac- 
counting for systematic uncertainties. A best fit sin 2 2a 
is determined for each Am 2 by minimizing the % 2 



function of sin 2a and these systematic coefficients. 



as a 
Fig- 



3 



ure g show the ratios of the measured rate of (is e ,17 e ) CC 
events to the Monte Carlo predictions minus 1. The in- 
ner errors are the statistical errors. In order to show the 
magnitude of the systematic errors, the outer errors in- 
clude all systematic errors added in quadrature (in the 
analysis all correlations are taken into account). The 
shaded area in the figure corresponds to the systematic 
error band from uncertainties in the predicted electron 
flux (primarily from the error in the branching ra- 
tio) and uncertainties in the measured flux at the detector 
(primarily from the r/3,772 difference, and the 2% error in 
the electron energy scale) . The curves correspond to the 
predictions for oscillations with sin 2 2a — 0.01, and Am 2 
of 100 and 1000 eV 2 . 
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FIG. 3. (a) Excluded region of sin 2 2a and Am 2 for 
— * v e oscillations from the NuTeV analysis at 90% con- 
fidence is the area to the right of the dark, solid curve, (b) 
NuTeV limits for V u —*V e . 



At all Ar 



the data are consistent with no observed 



are within one standard deviation of zero) . The frequen- 
tist approach []l3| is used to set a 90% confidence upper 
limit for each Am 2 . The limit in sin 2 2a corresponds to 
a shift of 1.64 units in \ 2 from the minimum (including 
all systematic errors). The 90% confidence upper limit is 
shown in Fig. |^(a) for — ► v e . Also shown are limits 
from BNL-E734 M and BNL-E776 pf. For sin 2 2a = 1, 



Am 2 > 2.4 eV 2 is excluded, and for Am 2 > 1000 eV 2 , 
sin 2 2a > 1.6 x 10~ 3 (the best fit is at 1000 eV 2 is 
sin 2 2a = (0.4 ± 0.9 x 10~ 3 ). In the large Am 2 region, 
NuTeV provides improved limits for — > v e oscillations. 

Similarly, the limit for — > ~D e is shown Fig. |3|(b). 
Also shown are the LSND || results and limits from 
KARMEN UJ. For the case of sin 2 2a = 1, Am 2 > 
2.6 eV 2 is excluded, and for Am 2 > 1000 eV 2 , sin 2 2a > 
1.1 x 10~ 3 (the best fit is at 1000 eV 2 is sin 2 2a = 
(— 0.3± 1.1 x 10~ 3 ). In the mode, our results exclude 



the high Am end of 



oscillations parameters 



favored by the LSND experiment, without the need to 
assume that the oscillation parameters for v and v are 
the same. These are the most stringent experimental 
limits fl6[ ] for v^ty^) —> v e (j> e ) oscillations in the large 
Am 2 region. 
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